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Abstract. The distribution of energy and momentum deposited by a fast parton in a medium 
of thermalized quarks, or the source term, is evaluated in perturbative thermal field theory. 
The calculation is performed by directly evaluating the thermal expectation value of the quark 
energy-momentum tensor. The fast parton is coupled to the medium by adding an interaction 
term to the Lagrangian. I show that this approach is very general and can be modified to 
consider more realistic modeling of fast parton propagation, such as a fast parton created in an 
initial hard interaction or the evolution of a parton shower due to medium induced radiation. 
For the scenario considered here, it is found that local excitations fall sharply as a function 
of the energy of the fast parton. These local excitations couple directly to the sound mode in 
hydrodynamics and are important for generating an observable Shockwave structure. This may 
have implications for the trigger pr dependence of measurements of azimuthal dihadron particle 
correlations in heavy-ion collisions. In particular, one would be less likely to observe a conical 
emission pattern for increasing trigger pr- 



1. Introduction 

The observation of the quark gluon plasma (QGP) created in heavy-ion collisions PQ is one of 
the most exciting discoveries of the past decade. Among the many interesting experimental 
results from the heavy- ion programs at the Relativistic Heavy Ion Collider (RHIC) and the 
Large Hadron Collider (LHC), the suppression in the production rates of energetic, or hard, 
leading particles and jets relative to a naive superposition of nucleon-nucleon collisions is one of 
the most striking [2JEH3]. These phenomena are often collectively referred to as "jet quenching" 
and are reflective of the interaction of hard partons with the QGP medium which they traverse. 

Another striking result from heavy-ion experiments is that the QGP formed at currently 
achievable center-of-mass energies may behave as a nearly ideal fluid [5]. In particular, the 
low energy, or soft particles created in heavy-ion collisions show a strong flow profile consistent 
with nearly ideal fluidity. The two results cited above lead to the natural question of how the 
QGP responds to a propagating fast (hard) parton. More generally stated, one would like to 
understand how the hard and soft physics of the QGP join together. 

In this proceedings I will attempt to make some progress in that direction. I will present 
a new approach for coupling the hard physics of fast parton propagation to the underlying 
thermal medium. The idea is to couple the fast parton to the QGP via an interaction term 
in the Lagrangian. The dynamics of the medium in the presence of the fast parton can then 
be directly obtained by taking the thermal expectation value of the energy-momentum tensor 



(EMT), denoted T^ u . Although I will here consider an asymptotically propagating fast parton, it 
should be stressed that the approach used in this proceedings can be easily extended to consider 
more realistic scenarios. For instance, one can modify the interaction term to consider back to 
back fast partons created in an initial hard interaction, which is more relevant to phenomenology. 
Or, equally as relevant, one can incorporate [6] the medium induced evolution of a primary fast 
parton into a full parton shower through radiative processes [7]. 

Ideally one would like to evaluate components of the EMT directly, but for the reasons 
discussed in [8] the focus here will instead be on the divergence of the EMT, d^T tiV . I will 
refer to this quantity as the "source term", in keeping with common terminology. The source 
term carries information about the rate of energy transfer to the medium and acts as a source 
(hence the name " source term" ) for the evolution of the underlying medium in the presence of a 
fast parton. By using an effective theory such as hydrodynamics, one can evaluate the medium 
response to a fast parton from the source term. A rigorously derived source term thus provides 
guidance for phenomenological studies of shockwave formation in the QGP |10j . 

This dual utility of the source term makes it a powerful quantity, particularly when applied 
to the realistic scenarios discussed above. For instance, the amount and distribution of energy a 
parton shower loses to the underlying medium can be obtained from the source term. This 
is useful in the description of jet observables [9], which are a more powerful probe of the 
QGP than leading particle observables. Additionally, such information may be useful in the 
experimental/jet medium background separation in heavy-ion experiments [3l 111], 

2. Formalism and Results 

The approach here is to begin with the EMT of the underlying medium. For simplicity I will 
consider a medium of massless quarks/antiquarks, leaving the inclusion of medium gluons for a 
forthcoming study. Thus, apart from the coupling strengths, the results will be the same as for 
a QED plasma. The EMT is given by [T2] 

= % -4> (y D" + 7 17 D^j V - 9^^, (!) 

where 

^ = ^0^,D^ = d^-igA^t a (2) 

and 

^^D u 7p = ^^D u ip-^j fl D* v ip. ( 3 ) 

In the above equations, g is the strong coupling, t a are the SU(3) generators in the fundamental 
representation, and conventional slashed notation is used, 4 = Tu^j etc. A summation over 
color, spin, and the active number of quark flavors is implied in the EMT. 

As mentioned in the Introduction, one must specify how the fast parton couples to the medium 
in order to investigate the medium response. The choice that is adopted here is to model the 
fast parton as an external color current which couples to the Lagrangian: 

se^se-A%%, (4) 

where for the moment I do not specify the explicit form of j (j u should not be confused with 
the source term, J v ). Non-Abelian gauge symmetry is preserved by the replacement made in 
© if D?j£ = 0. 



Figure 1. The Feynman diagrams contributing to (d^T^) in the presence of the interaction 
term, j%. The contribution from the energy- momentum tensor and external current insertions 
are discussed in the text (see also Figure [2]). The diagram in Figure 1(a) can be traced back 
to terms in the energy- momentum tensor (see equation ([1])) which go as ^dip, whereas the 
diagram in Figure 1(b) originates from terms that go as gip'yA'ip. 



In the presence of the interaction term, ja , the lowest order Feynman diagrams needed to 
calculate the thermal expectation [d^T^) are shown in Figure [TJ The diagram in Figure 1(a) 
arises from terms in ([T|) which go as tp^dtp, whereas Figure 1(b) arises from terms that go as 
gtp^Aijj. Note that a two gluon exchange is necessary to couple to the EMT, which is a color 
singlet quantity (a two photon exchange is also necessary in QED, from Furry's theorem). The 
convention used here is that the standard Feynman quark-gluon vertex contributes ig^t a and 
the Feynman gauge is used for gluon propagators. 

One must also determine what the correct Feynman rules for the EMT and external current 
are. To determine these non-standard contributions, which are isolated in Figure [21 one first 
assigns the appropriate momentum to each derivative. Recalling that one is here interested in 
{dnT^ u ) and that the final result is Fourier transformed into position space, it is found that the 
value of Figure [2(a) is [8] 

■ p -ix-(pi-p 2 ) , x 

4 ((p! ~pI)Y +J*(3jf 2 +jg -P^f, +t 2 )) ■ (5) 

One applies the same procedure to Figure [2(b) where I choose the convention that the gluon 
momentum flows away from the external current (or into x and any vertex) in all cases. The 
result is [8] 

_ i9e -., (w+K - K , to + P3 _ (6) 

which is only valid when the gluon field in Figure [2(b) connects with the source. 
Finally, for the external color current in Figure [2(c) , one has very generally 

-i f d*zj«(z)e iz ^ (7) 

for the case of an external current which contains one power of g. As mentioned in the 
Introduction, the focus here will be on an asymptotically propagating fast parton. In that 
case ja = gQ a (t)U fl 5 3 (z — ut) where u is the fast parton's velocity and = (l,u). gQ a (t) 
is the charge of a classical particle in QCD (defined by Qf Qj = 5ijC2i, with C21 the quadratic 
Casimir in representation % (3 for a gluon, 4/3 for a quark)) and Q a (t) evolves in time according 
to Wong's equations [32]. The time dependence of Q a (t) is of higher order in g and will not be 
considered here. For this choice of ja Figure [2(b) simplifies to 



-2iri g Q a U a 5(p A • U). 



(8) 




Figure 2. Contributions to the Feynman diagrams of Figure [T] which come from inserting the 
energy-momentum tensor and external current interaction term. Their values and how they are 
obtained are discussed in the text, see specifically equations © - (J7J). 



As discussed in the Introduction, the approach used in this paper can easily be extended to 
consider more realistic scenarios than the single asymptotic parton considered here. One can 
modify the interaction term to consider back to back fast partons created in an initial hard 
interaction or incorporate the medium induced evolution of a primary fast parton into a full 
parton shower through radiative processes. 

To obtain the Green's function structure of Figure[T]I use the conventions of real time thermal 
field theory as outlined by Das [13] • The details of the calculation will not be presented in this 
proceedings, but will be presented in a follow-up publication in which medium gluons will also 
be considered. Using standard Feynman rules for real time thermal field theory in addition to 
the special rules presented in equations dSJ) - ([TJ) the result for the source term is obtained as [8] 

(d,T^{x)) = -4iN F C 2 g 4 J d ' n ^f n e-^ + ^5(p 2 3 )n F (p 3 )G R (p 4 )G R ( P3 + p 4 )G R (p 5 ) 
x S(p 4 ■ U)5( P5 ■ U) [2(p 3 • Ufvl ~ U 2 P3-P4P U 5 - U v (P3 ■ U)(2p 3 -p 5 +p 4 -p 5 )] 

(9) 

where G R (p) = (p 2 + iep )" 1 is the retarded Green's function, tif(p) = {e [P °^ T + l)" 1 is the 
Fermi distribution function, T is the temperature, and Np is the number of active flavors. 

Equation (jUJ) is the central result of this proceedings. One can obtain an intuitive picture 
of the underlying physics it contains by considering the Green's function structure. Reading 
from left to right in the first line it is seen that a particle from the heat bath (represented by 
the thermal distribution 5{p\)np{p 3 )) absorbs a gluon from the external current (represented 
by Gr(p4)) and then continues to propagate (shown by G R {p 3 + Pa)) until it absorbs another 
gluon from the external current (represented by G R {p^)). A final propagator representing the 
reabsorption of the particle by the heat bath (which would be given by G R (ps +P4+P5)) has 
been cancelled by the momentum structure of {d^T^"). 

It is useful to consider certain approximations and limits of equation Q. A kinetic theory 
calculation of the source term generated by an asymptotically propagating fast parton was 
performed in |15j . Within the hard thermal loop (HTL) approximation equation (|9j) should 
reduce to the kinetic theory result (the HTL approximation is formally equivalent to the Vlasov 
equation [E]). The interested reader can verify (as this author has done) that (|9|) indeed 
reproduces the result of |15| . This verification gives confidence to the approach introduced in 
this proceedings. 

I now consider an application of ([9]) by considering the source term ansatz discussed in |17j . 
The authors of [T7] proposed the following simple form: 

jp 

(d^ix)) ~-jr{U v - \d v ) J(x - ut). (10) 



The utility of ()10p is that one can encode basic features of the source term in a compact way. 
The dimensionfull coefficient A parameterizes local contributions from the source, that is, terms 
which globally integrate to zero (and for instance do not contribute to dE/dt) but may still be 
important for exciting the medium. It was found in [T7] that within linearized hydrodynamics 
a double peaked structure in the azimuthal emission spectrum associated with the source given 
by (|10p only appeared for rather large values of A (on the order of 0.5 fm or higher for 20 GeV 
total energy deposited into the medium). In what follows, I will extract A from equation Q. 
To extract A, one needs the relation [5] 

1 r* 



X = m d 3 xx{d^ x (x)) (11) 

dt J 

where I have chosen to use the x component of the source to obtain A (according to the ansatz 
of (jlOp any component would work). dE/dt here is the collisional energy loss rate, which can be 
obtained from the source term as 



dE 
~dt 



d 3 x(d fl T> l0 (x)). (12) 



For this proceedings I will use the HTL limit for the collisional energy loss: 

dE _ m 2 D C 2 a s A tanh" 1 [u] \ ^ p max 
dt htl 2 V 7 2 it J p min 

where 7 = (1 — u 2 ) -1 / 2 and based on the HTL approximation and physical reasoning p ma x ~ T 
and Pmin ~ triD- For simplicity, I will use the relativistic limit (7 3> 1) of (jT3j) . 

With dE/dt in hand, A can be obtained directly from ([9]) by using (jlip and an integration 
by parts of the general form: 

j d 3 xd 3 pxe ip - x f(p) = i J d 3 xd 3 pe ip - x d p J(p). (14) 

Care must be taken in using (|14p because A has an infrared divergence. The simplest way 
around this, and what I will adopt here, is to introduce a mass term in the Green's functions 
Gr{pa)Gr{p^) in ([9|) given by the thermal gluon mass m 2 = m 2 D /2. 

It is also necessary to consider carefully the momentum scales involved in the integration. 
For consistency, one should only apply the full result of equation ([9]) down to some momentum 
scale \pi\ ~ q* where gT C 9* < T. When \p^\ ~ gT the full result contains contributions from 
higher order terms and only the HTL approximation is consistent. From a practical point of 
view in the calculation of A this means dividing terms in the integration such that 

A = A F when |p 4 | > q*, 
A = Ahtl when \p$\ < q*, 

where Af is obtained from ([9]) and Ahtl is from taking the HTL approximation of ([9|). The 
final result should be independent of q* (this type of analysis for the separation of scales in 
calculations of finite temperature field theory was introduced in |18j). 

With these technicalities in place, I will present the result for A and further show that it 
is independent of q*. The result is presented in Figure [3] as a function of 7 = (1 — u 2 )~ 1 ^ 2 
for 5 = 1 and T = 400 MeV. It is interesting that A falls rather sharply as a function of 7. 
As was mentioned above, A parameterizes local contributions from the source and is important 
for exciting the medium. In particular, A couples directly to sound modes (and not to diffusive 
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Figure 3. In the source term of equation (|10p A couples directly to sound modes and is 
important for generating observable Mach cone signals. The contribution to A from the source 
term ([9]) is plotted above for g = 1 and T = 400 MeV [8]. It is clear that A drops as a function 
of 7 = E/m which could have implications for the trigger px dependence of azimuthal dihadron 
correlation measurements (see discussion in text). The inset shows that the result for A is largely 
independent of the separation parameter, q* (see text for details). 

modes) when (|10p is used as a source term for hydrodynamics. It was found that large values of A 
were crucial to the appearance a conical Mach-like emission spectrum in |17j . The experimental 
implication of the dependence of A on 7 as shown in Figure [3] could be found in the trigger px 
dependence of measurements of azimuthal dihadron particle correlations. In particular, a conical 
emission pattern would be less likely to be observed for increasing trigger px, which indeed seems 
to be the case [19] . The inset of Figure [3J shows that the result is largely independent of the 
separation parameter, q*, introduced above. 

3. Conclusion 

In this proceedings I have presented a new approach for coupling the hard physics of fast parton 
propagation to the underlying thermal medium. The principle idea is to couple the fast parton 
to the QGP via an interaction term in the Lagrangian. The dynamics of the medium in the 
presence of the fast parton are then directly obtained by taking the thermal expectation value 
of the energy-momentum tensor. 

Using this approach the source term for a medium of thermalized quarks in the presence of 
an asymptotically propagating fast parton was obtained and presented in equation ([9]). It was 
then found that local excitations, which may be important for generating observable Mach cone 
signatures, fall sharply as a function of 7, as shown in Figure [31 This may have implications 
for the trigger px dependence of measurements of azimuthal dihadron particle correlations in 
heavy-ion collisions. In particular, a conical emission pattern would be less likely to be observed 
for increasing trigger px, which may indeed be the case [19j . 

A future publication will include medium gluons and present the details of the calculation 
|20j . The approach introduced here is general and can be modified to consider the medium 
response to tagged jets [21] or medium induced full parton showering constrained by realistic 
radiative energy loss calculations |22j . 
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